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Abstract: We claim that factorization implies that the evolution kernel, defined by 
the logarithmic derivative of the N-th moment of the structure function dlnF^ /dlnQ 2 , 
receives logarithmically enhanced contributions (Sudakov logs) from a single source, 
namely the constrained invariant mass of the jet. Available results from fixed-order cal- 
culations facilitate Sudakov resummation up to the next-to-next-to-leading logarithmic 
accuracy. We use additional all-order information on the physical kernel from the large- 
(3q limit to model the behaviour of further subleading logs and explore the uncertainty 
in extracting a s and in determining the magnitude of higher-twist contributions from a 
comparison with data on high moments. 



1 Introduction 



The standard collinear factorization in deep inelastic scattering (DIS) at twist two is 
based on separating hard physics associated with scales of order of the momentum trans- 
fer Q 2 into coefficient functions, and softer physics associated with the structure of the 
target into operator matrix elements which have the interpretation of parton distribution 
functions. 

Considering the kinematic limit x — ► 1 (where Bjorken x is defined by x = —q 2 /2pq = 
Q 2 /2pq), an additional, well-separated scale emerges, W 2 = Q 2 {l—x)/x which represents 
the invariant mass of the hadronic system. This is the characteristic scale for dynamics 
of the jet in the final state. When W 2 is much smaller than Q 2 , perturbative corrections 
(Sudakov double logs) related to soft and collinear radiation which form the jet dominate 
the coefficient functions. These corrections become so large that the perturbative ex- 
pansion breaks down. Consequently, they must be resummed to all orders. In standard 
analysis of structure function data the difficult large- x region is usually avoided simply 
by putting a lower cut on W 2 . For example, in the MRST analysis |1] a cut is put such 
that W 2 > 12.5 GeV 2 . 

To describe the structure functions at large Bjorken x, it becomes natural to further 
factorize the process and treat the jet separately. Such factorization can indeed be 
accomplished @-|§, facilitating the resummation of the logarithmic corrections. A 
schematic picture of large- x factorization is shown in Fig. [I| 




Parton distribution 



Figure 1: Factorization of DIS structure functions at large Bjorken x. 



It is well known that higher-twist corrections are also enhanced at large x: they 
appear as powers of A 2 /W 2 rather than A 2 /Q 2 . At sufficiently small W 2 the Operator 
Product Expansion (OPE) itself breaks down and needs to be resummed. Given the 
complexity of this expansion P|-||12||, the only hope this could ever be achieved is if 



a small subset of matrix elements dominates at each twist, simplifying greatly the ex- 
pansion. There are two ways in which the leading twist simplifies at large x. Firstly, 
the valence quarks dominate over the gluons and the sea. Secondly, as explained above 
(Fig. |l|) the dynamics of the outgoing jet, associated with the scale W 2 , decouples from 
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the rest of the process. The question then arises whether any of this holds beyond the 
leading twist. 

Recently, a first theoretical analysis of the large- x limit of twist four has been per- 
formed fl3| , suggesting that indeed, this kind of simplification occurs to all orders in the 
twist expansion, and thus its resummation at large x may be possible. 

Ref. JHJ first establishes, in the concrete case of F 2 , the cancellation of the infrared 
renormalon ambiguity at twist two with the quadratically divergent contribution of twist- 
four operators, associated with their mixing with twist two. Renormalons reflect the 
ambiguity in separating contributions of different twists. Independently of how this 
separation is implemented, within twist four there is a twist-two like ingredient, which 
is proportional to the twist two matrix elements. The question that arises is what is 
the significance of this twist-two like component compared with other "genuine" higher- 
twist contributions. Renormalon-based models for power corrections (14[]-[18| assume its 
dominance whereas other approaches may neglect it altogether. 

The conjecture of |13[ is that at large x the twist-two like ingredient within the higher 
twist indeed dominates. Analysis of the leading-order coefficient functions of twist four 
for i*2 and Fl provides a physical picture which is consistent with this conjecture. The 
final state which dominates the twist-four contribution at large x is the same as the 
twist-two final state: in the case of F 2 this state contains a single energetic quark which 
initiates the jet. Moreover, the analysis of also reveals that the multi-parton initial 
states at twist four, which distinguish it from twist two, are dominated at large x by 
particular configurations which are twist-two like: a single quark in the initial state car- 
ries most of the momentum, whereas additional gluons carry small momentum fractions. 
Consequently, the corresponding twist-four matrix elements effectively depend only on 
the light-cone separation between the quark fields, similarly to twist two. 

The conclusion is that factorization, similar to the one described by Fig. [1], may hold 
beyond the perturbative level. The essential difference is that beyond twist two, gluons 
of virtuality of the order W 2 are exchanged between the jet and the remnants of the 
target. Their effect, to all orders in the OPE, can be taken into account through a 
single "shape function" of N/Q 2 , which multiplies the twist-two contribution in moment 
space. This way, the conjecture of ref. [jnj translates into a concrete understanding of 
how higher-twist effects should be parametrized at large x. 

The need to resum power-corrections to all orders as a consequence of kinematic 
thresholds exists in a wide range of hard processes, including, for example, Drell-Yan 
pair production and event-shape distributions. In spite of the different nature of these 
processes, the analogy is quite useful. In particular, the appearance of a new non- 
perturbative distribution |19[, the shape function, near the kinematic threshold, is com- 
mon to all these processes, and so is the relation between the resummed perturbative 



distribution and the properties of this non-perturbative function [20, 21, B2|. 



It should be stressed that the proposed non-perturbative factorization formula for DIS 
structure functions has not been derived from the OPE. It requires that the Sudakov 
resummed coefficient function of higher-twist contributions to F 2 would be the same as 
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that of twist two. Moreover, the large- N behaviour of higher- twist anomalous dimensions 
should coincide, asymptotically, with that of twist two. This highly non-trivial structure 
is yet to be verified by explicit calculations. 

Having assumed the dominance of these higher-twist contributions which mix under 
renormalization with the leading twist (and are therefore proportional to the latter), 
renormalon resummation becomes absolutely essential for any higher-twist analysis. In- 
deed, the separation between the leading and higher twist is ambiguous. But, more 
importantly, prior to introducing any parametrization of power suppressed contributions 
one must be sure that any large perturbative corrections have already been taken into 



account |23|, 21| . This concerns, in particular, running-coupling (or renormalon- related) 
perturbative corrections which are always parametrically larger than the corresponding 
ambiguity, and thus, by the previous assumption larger than the higher-twist contribu- 
tion itself. Failing to take these perturbative corrections into account, the parametriza- 
tion of the higher twist becomes meaningless. For example, the values of the extracted 
"higher twist" parameters will strongly depend on the order of truncation of the pertur- 
bative expansion, the renormalization scale used, etc. 

The unresolved issue of higher twist in DIS structure functions dates far back. It was 
understood very early on that higher-twist effects are especially important at large x. The 
limitations have always been both insufficient data and lack of theoretical understanding 
of this region. In the recent years there has been much activity in this field, primarily 
owing to progress in perturbative calculations and resummation. Several groups have 
performed fits to data p^|-[^|, incorporating (or not) soft gluon resummation and 
parametrizing the higher twist in various ways. What lacks however, is a physical picture 
that stands behind these parametrizations. In particular, one is bound to ask eventually 
how the higher-twist parameters are related to operator matrix elements within the 
OPE, namely which multi-parton correlations are being measured. Ref. |13j suggests an 
answer[} More theoretical and experimental work is required to establish it. Here we 
make a further step in this direction. We concentrate mainly on perturbative aspects, 
constraining further the structure of the Sudakov exponent and examining where we 
stand on the long track to power accuracy. 

As mentioned above, large-x kinematics singles out very specific radiative correc- 
tions which dominate the coefficient function: these are the logarithmically-enhanced 
terms associated with the evolution of the jet. The jet function can be written as an 
exponential in moment space, where the exponentiation kernel at a given order can be 
systematically deduced from the standard loop expansion of the same order. To a first 
approximation (leading-log accuracy) the kernel simply corresponds to a single gluon 
emission calculation, whereas multiple emission is accounted for by the exponentiation 
(or solving the evolution equation with that kernel). Refining the calculation of the ker- 
nel to next-to-leading logarithmic accuracy amounts to taking into account the leading 
contribution associated with the running of the coupling. A natural procedure to fur- 

1 A previous suggestion with some similarity is that of ref. fl33[. 
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ther improve the approximation to the kernel is a renormalon calculation with a single 



dressed gluon [^XJ, £!| . In this way some running-coupling effects are taken into account to 
any logarithmic accuracy, making the exponent renormalization-scale invariant. Conse- 
quently the related power-corrections can be systematically parametrized, in accordance 
with |3| . 

The essential difference between dressed gluon exponentiation (DGE), and other ap- 
proaches to Sudakov resummation, is that it takes into account some all-order informa- 
tion on the exponentiation kernel itself. However, since the renormalon calculation is 
restricted to the large- Nf (or, equivalently, large (3q) limit, it is necessary to combine it 
with what is known about the first few orders in the kernel from fixed-order calculations. 
One of the goals of this paper is to combine the large-/?o result |22j for the Sudakov 



exponent of F2 with the state-of-the-art knowledge of the MS anomalous dimension and 
coefficient function. The latter practically^ allows a complete |33] next-to- next-to-leading 
logarithmic (NNLL) accuracy calculation, so together with the all-order large- Nf result, 
it seems that the Sudakov exponent as a whole is well constrained. Nevertheless, we will 
see that without making further assumptions concerning the structure of the Sudakov 
exponent, there still is a significant uncertainty concerning the magnitude of sub leading 
logs. This has important consequences for phenomenology. 

In this paper we apply, for the first time, DGE with the corresponding parametriza- 
tion of higher-twist corrections to DIS data. We restrict our attention here to the large-x 
region, by analysing moments N > 5. The moment space analysis is advantageous from 
a theoretical point of view in several respects. First of all, it facilitates the implemen- 
tation of target-mass corrections through the use of Nachtmann weights. In addition, it 
simplifies significantly the resummation formulae for the coefficient function as well as 
the higher-twist corrections. 

The paper is organized as follows. The theoretical part (section 2) begins by recalling 



the non-perturbative factorization formula of [13|. We then return to discuss factoriza- 



tion on the perturbative level and its consequences for resummation (section 2.1). We 
emphasize in particular the fact that Sudakov logs in the physical evolution kernel are 
exclusively related to the jet function and not to the soft function. This translates into a 
more predictive formulation of Sudakov resummation than that commonly presented in 
the literature. This discussion is followed by an all-order analysis of twist-two at large N, 
first (section 2.2) for the physical kernel, and then (section 2.3) for the coefficient func- 
tion. Sections 2.4 and 2.5 are devoted to combine the all-order large- Nf information with 
the available coefficients in QCD and discuss in some detail the uncertainty involved in 
such a procedure. Section 3 is devoted to data analysis. Matching the DGE exponent 
into the known |55[|- P0[] next-to-next-to-leading order (NNLO) result we obtain an im- 
proved perturbative prediction for the scaling violation at large N, which we use as a 
baseline for the study of power corrections. The experimental moments N = 4 - 11 of F 2 



2 The non-Abelian contribution to the MS anomalous dimension is still not available analytically as 
a function of N. However, there is a reliable estimate j}8| of the leading In N contribution based on the 
computation of specific moments |37j . 
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for the proton are calculated based on SLAC and BCDMS data. Each of the moments 
is fitted separately, to extract the valence quark distribution, a s and the magnitude of 
the higher twist. We use the stability of the extracted value of a s as a function of iV as 
a measure of the quality of our theoretical description of F%, and compare the results to 
the standard NNLO analysis. In section 4 we summarize our conclusions. 



2 Factorization and the Sudakov exponent by DGE 
2.1 Factorization 

The non-perturbative factorization formula [13] for F2 at large N takes the form: 

F 2 N (Q 2 ) = £ dxx N - 2 F 2 (x, Q 2 ) = H (Q 2 ) J n (Q 2 ; /xf.) q N (fi%) J NP (NA 2 /Q 2 ) . (1) 

We will assume that it holds up to corrections of order 1/N. Here, the first three 
factors correspond to the standard, perturbative large- 2 factorization [@]-0: H (Q 2 ) 
is the hard part of the coefficient function depending on the momentum transfer Q 2 , 
Jn (Q 2 ;fip) is the Sudakov-resumed jet function depending primarily on the invariant 
mass of the jet Q 2 /N, and g./v(/-40 is the twist-two quark matrix element. The last factor 
jnp ^jVA 2 /Q 2 ) = 1 + Ki + • • •, resums the dominant higher twist corrections to all 
orders [TB[ . 

Some important clarifications are due concerning the perturbativeQ factorization for- 
mula. First, the coefficients in the hard function H (Q 2 ) are finite at large N, so this 
function does not play any role in the large- N limit. The jet and the quark distribution 
are defined as follows: 

■/(Q 2 (l-0/£;/4) = ^Im/ d'ye-^^y (0|Tr(*(£)*(0)^)|0), (2) 

with £ ee x/z, and 

q(z;i&)= [°° ^e^Wp *(y)^(0) p) (y 2 = 0) (3) 

v ' J-00 py x 1 

where the gauge in (Q) and ([|) is chosen as A_ = and A + = 0, respectively. Here 
the lightcone '+' direction is defined by the incoming hadron p, and the lightcone '— ' 
direction by the outgoing jet, (q + xp). J is simply a projection of the quark propagator 
in the axial gauge pA = 0. At leading order \l/(y)\I / (0) is the free quark propagator and 
J = 8(z — x); at higher orders J depends on (q + zp) 2 = Q 2 (l — — Q 2 (l — £) but 
also, through the gauge fixing, on the hard scale pq. The latter dependence is unrelated 
to logarithmically enhanced terms and it can be neglected. Corrections on this scale 
appear through the hard function H{Q 2 ). 

3 The reader in referred to 03] for discussion on power corrections. 
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In spite of the explicit dependence on the gauge, J (Q 2 (l — O/CjA'f) an d Q.{ z \1*f) 
are gauge invariant. In case of different gauge choices, an appropriate path-ordered 
exponential should be introduced. For J the path goes along the incoming quark di- 
rection (p), connecting the point y to infinity and then infinity to 0, whereas for q 
the path goes along the jet direction between y and 0. The different paths are shown 
by the dashed lines in figure 0. The functions in (|IJ) are defined in moment space, 
Jn(Q 2 ', A*f) = Jo dx x N ~ l J ((1 — x)Q 2 ; n 2 F ) and <7jv(/x|0 = Jq dx x N ~ x q(x; n 2 F ), where 
convolution in the momentum fraction simply becomes a product. 

y 

j ; : 
q 

o y 



V ! i 

Figure 2: Factorization of F 2 at large x. Full and dashed lines stand for a dynamical 
quark and a Wilson line, respectively. The vertical lines correspond to the incoming 
lightcone direction p and the horizontal line to the outgoing jet. The three pictures 
correspond to the jet, the quark distribution and the soft function, respectively. 

The fj/p dependence in (Q) and (|) is there to remind us that the definition of these 
matrix elements requires a prescription (factorization scheme and scale) to deal with 
infrared and ultraviolet divergence, respectively. Since F 2 is a physical quantity, this 
dependence must cancel out in the product. 

The standard formulation of large-x factorization ||-[f| includes an additional soft 
function V: 

F 2 (x, Q 2 ) = H (Q 2 ) J ((I - x)Q 2 ) V (x- /4) ® q{x; /4). (4) 

defined 0, ^ [IT], |] by the path-ordered exponential shown in the lower picture in fig- 
ure 0, which represents the DIS process in the Eikonal approximation. The proof of 
factorization relies on separating a generic Feynman diagram into subprocesses, each 
depending on a single scale. In inclusive DIS there are two external scales: Q 2 and W 2 = 
Q 2 (l — x)/x ~ Q 2 /N. The corresponding subprocesses are H{Q 2 ) and J{Q 2 {\ — x)). 
Momentum scales which are parametrically smaller than W 2 , e.g. 0(Q 2 (1 — x) 2 ), are 
accounted for by the soft subprocess V(x, p, F ). 

It is important to note that V(x, fi 2 F ) does not depend^ on any external scale in 

4 Dependence of this function on Q 2 or on Q 2 /N 2 is sometimes introduced through the factorization 
scale. This, however, can be misleading, and we will avoid it in this formal discussion. 
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the DIS process, and it is therefore possible and useful(!) to eliminate it from the 
factorization formula Although not written explicitly in eq. in order to define 
the soft function one actually performs two factorization procedures with two different 
scales fi Fl 2 : one in the definition of J ((1 — x)Q 2 ; fi Fl ) and the other in the definition 
of q(x; n F2 ). This can be intuitively undersood from Fig. [I], where the soft blob is attached 
to jet and to the quark distribution. In any case some component of the soft blob is 
included in the jet and some in the quark distribution. Since V represents the evolution 
of the quark distribution as well as that of the jet function with /j, F , it can be eliminated 
altogether by choosing the same scale in both. This is demonstrated in figure |||. The 





Quark distribution 



Figure 3: Factorization of F2 at large x. The soft function is split between the jet and 
the quark distribution. Different separations amount to evolution of q(x,fi F ). In our 
formulation the upper 'soft' blob is understood as part of the jet while the lower one as 
part of the quark distribution function. 

factorization formula is then 

F 2 N (Q 2 ) = H (Q 2 ) J n (Q 2 ;/^) V n (/4,/4 2 ) q N (»%) 

= H(Q 2 ) J N (Q 2 ;v 2 F2 ) q N {n 2 F2 ). (5) 

where 

Jn (Q 2 ; /4 2 ) = j n (Q 2 ; /4J v n (/4u a4 2 ) • 

In the factorization formula (H) the soft function appears only implicitly through 
the dependence of Jn(Q 2 ',^f) and qN^p) on n 2 F . This function has some remarkable 
properties which will become relevant in the following. As discussed in ref. [|] (see 



section 7 there) and in |TI|, |5| its exponent contains just a single In N at any order in the 
coupling. It thus follows from the factorization formula that apart from these /z^-related 
single logs, the only source of logarithmically enhanced contributions to the F2 coefficient 
function is in the constraint on the invariant mass of the jet. This is true to all orders in 
perturbation theory. As we discuss further in the next section, this implies, in particular, 
that the jet function J((l — x)Q 2 ; fi 2 F ) alone determines the Sudakov exponent in the 



physical kernel ciln Ff \Q 2 )/d In Q 2 
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We stress that the absence of additional sources of Sudakov logs is strictly related 
to the inclusive nature of structure functions. Note that, due to different kinematic 
circumstances, sensitivity to large-angle soft emission which is associated with the scale 
Q 2 /N 2 does generate Sudakov double logs in other observables. Examples are provided 
by Drell-Yan pair production and event-shape distributions. 

Finally, in the notation of ||, |39|, |33| , our statement that the only source of Sudakov 
logs in the physical kernel is in the jet function means that D^ ls vanish to all orders. 
Indeed, recently, it was shown [|4]], based on explicit calculations in MS that -DJ? IS = 



D® ls = 0, although the conviction that D Dls vanishes to all orders was not made. 



2.2 The physical kernel 

It is instructive to formulate first Sudakov resummation in the physical kernel defined by 
the logarithmic derivative din (Q 2 ) / dlnQ 2 . We claim that the latter can be written, 
to any order in perturbation theory, as 

d\nF 2 N (Q 2 ) d\nJ N (Q 2 ^ 2 F ) C F r l x N ^ -1 , , 
dlnQ 2 dlnQ 2 Po Jo 1-x V 1 XN 

= 9JL rduT(u)(Q 2 /A 2 Y U (N u -l)T(-u)BlJC}(u), (6) 
Po Jo v ' 

where, as before, terms that are finite at large N are neglected. Here we used the scheme 
invariant Borel representation |[42||, 



roc , . _ u 

/C(^ 2 ) = ^ duB[K](u)T(u) (/i 2 M 2 ) , (7) 
and T(u) is the Laplace transform of the coupling^], 

A(Q 2 )= / duT(u)(Q 2 /A 2 ) . (8) 



In the following we will use the 't Hooft coupling A = a s (Q 2 )j3o/ii, defined by 

d A 

-A 2 (1 + 5 A); 5 = A//3o 2 , (9) 



so 



The function 



dlnQ 2 

T(u) = (u5) uS exp(-u5)/r(l + u5). (10) 



2 3 
r It ~ U ~ U 



B[K](u) = 1 + h- + k 2 - + k 3 - 



5 Eq. (^) is written in a renormalization scheme invariant way. Both T(u) and the coefficients of 
B[K,]{u) would change depending on the renormalization group equation for the coupling, such that the 
product remains invariant. 
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is the Borel representation of an N -independent effective charge, which is a function of 
the coupling having the expansion 



K = A + k x A 2 + k 2 A 3 + 



Our notation here is adopted for the large-/?o limit, where 5 = 0, T(u) = 1 and k n = k n . 
Beyond this limit the general relation can be written as an expansion in S, 

~k n = k n - 5 k n _ x n (*(n + 1) + lE - 1) + 0{5 2 ). (11) 

Like the physical kernel itself, /C is does not depend on any arbitrary factorization 
procedure. Note that even if K, is free of any infrared renormalons (as is the case in 
the large- (3 limit, and perhaps also in the full theory - see below), these do appear 
in d In F^ (Q 2 )/d In Q 2 through the factor (N u — 1) F(—u) which is associated with soft 
and collinear radiation under the large-x phase-space constraint p2|. This factor is 



understood as follows: gluon splitting off a quark gives rise to a 1/(1 — x) + singularity. 
The phase-space constraint determining the maximal gluon virtuality is k 2 < W 2 = 
Q 2 (l — x)/x. Using the Borel representation of the dressed gluon propagator, k 2 is raised 
to a power —u. Performing the phase-space integration, phase-space constraint promotes 
the x — ► 1 singularity to a cut: (1 — x)~ l ~ u . Going to moment space and neglecting 
0(1/ N) contributions one obtains: (N u — l)T(—u). As a result, at any positive integer 
value of u where B[K](u) does not vanish, the physical kernel will have a renormalon 
singularity. 

The all-order B[K](u) in the large- Nf limit has been computed in P2|, and it is given 

by 



B[K](u) 



1 sin7ru 



jCu 



large j3 



71U 



U 



+ 



1 



u/2 



(12) 



where C = 5/3 is the constant from the MS renormalization of the fermion loop0. 
Additionally, the first two coefficients k± and k 2 are known (see e.g. |M EO, 131]), by 



combining results for the coefficient function and the anomalous dimension in the MS 
factorization scheme. The relations will be given below. 

Examining the (large Nf) renormalon structure of the physical kernel @ with (|T2| ) we 
deduce that non-perturbative corrections at large iV appear as [^2[ [13| a multiplicative 
factor in moment space, as summarized by eq. ([l]), with 



J NP (NA /Q ) = exp { -Ui^f- 



C F ( NA 2 



Po V Q 2 



1 C F f NA 

2 U2 ^ Vo 2 " 




(13) 



This power-correction model effectively resums the twist expansion, under the assump- 
tion that the higher-twist contributions which dominate at large N, are those which 
mix with the leading twist. Fixing the parameters iO\ and uo 2 requires the knowledge 



6 We assume that A in 
equation is M. 



is defined in this renormalization scheme, while the renormalization-group 
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of the dominant contribution to the matrix elements of the corresponding higher-twist 
operators |13| . Clearly, the non-perturbative contribution J NP (NA 2 /Q 2 ) has a meaning 
only within a given regularization prescription for the renormalons in the perturbative 
jet function (Q 2 ; fi 2 F ). 

In the large-/?o limit /C is free of any renormalon ambiguity. B[JC](u) is an analytic 
function at any finite u, so renormalon singularities in the physical kernel are exclusively 
related to the end-point singularity at x — > 1. If K, is free of infrared renormalons 
also in the full theory, it follows^ that at large N the anomalous dimension of higher- 
twist operators (which mix with twist two), coincides with that of the twist two, as 
was already anticipated! in f!3|| . Another possibility not yet excluded is that beyond 
the large-/?o limit, B[K](u) would have infrared renormalons. For example, the poles 
at u — 1 and u = 2, which are exactly compensated by the sine factor in flT2"|), might 
become branch points. Such a cut structure would imply that the large N asymptotics 
of higher- twist anomalous dimensions differ from that of the leading twist, and a more 
complicated pattern of power corrections at large N would emerge. 

Finally, let us return to the last point we made on section 2.1. The absence an addi- 
tional source of Sudakov logs of the form D DIS ((1 — x) 2 Q 2 ) in inclusive DIS (such a term 
was included, for example, in can also be understood from general consideration 
on power corrections. In DIS the external Lorentz invariants are Q 2 and W 2 ~ Q 2 /N. 
On these grounds alone one expects power corrections to appear on either of these two 
scales and not on any other scale. For example, corrections such as A 2 N 2 /Q 2 should 
not appear. On the other hand, had there been a function P DIS ((1 — x) 2 Q 2 ) in eq. (|6|), 
renormalons associated with the scale Q 2 /N 2 would have appeared^ upon performing 
the integration over x, similarly to the situation in Drell-Yan production [ pE9"| , |2"2fl . As 
we argued in section 2.1, the perturbative factorization of structure functions does not 
require to introduce such a function and this problem is avoided altogether. 



2.3 The jet function 

Going from the physical kernel to \nF^ itself, the differential equation @ should be 
supplemented by an initial condition. Within the OPE, the initial condition is defined as 
the matrix element of the twist-two quark operator with N covariant derivatives, Qn{^ 2 f)- 
The operator has a non-trivial renormalization, which takes the form 

d^N(^) = _Cr lnNA 
a In \ip p 



7 As usual p3l 0], the nature of infrared renormalon singularities is related with the anomalous 
dimension of higher-twist operators. 

8 Let us recall that the anomalous dimension of twist three is known to coincide asymptotically with 
that of twist two |45| . 

9 The corresponding Borel integrand includes the factor (iV 2u — l)r(— 2m), so unless B[D](u) vanishes 
at all integer and half integer values of u, there would be power corrections on the scale Q 2 /N 2 . 
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up to corrections which are finite at large N. The effective charge 



A = A + a 2 A 2 + a 3 A 3 + ... 

is associated with the soft function V in section 2.1. It corresponds to the 1/(1 — x) + 
singular part in the Altarelli-Parisi splitting function, which can be computed using the 
Eikonal approximation, or the Wilson lines shown in the lower picture in Fig. It has 
been noted that A is a rather universal object HI], 0]: it is the anomalous dimension of 
a Wilson line with a cusp. Upon integrating @ and (|14]) and using the factorization 
formula ([!]) to write In Jjy as the difference between In and In (both regularized 
by the Borel variable m ^ 0), we obtain the following jet function: 



In J N (Q 2 ;n 2 F ) = ~ r duT(u) (Q 2 /A 

v ' On JO v 



r( _) ( iv- - 1) + rMtrrsm*,, 

11. V / 11. 



U K ' u 

where, as before we use the scheme invariant Borel representation, where 



(15) 



2 

Ii Ii 

B[A](u) = I + ~a 2 - + ~az— + . . . . (16) 

Note that the second term in the squared brackets cancels the u = singularity of the 
first, so in In Jjv (Q 2 ] fJ/p) the integral over u exists. This is in contrast with lng^r and 
lnF^, which are not perturbative entities. The dependence of \n J N (Q 2 ; fi 2 F ) on the 
factorization scale \xp and on the factorization scheme is guaranteed to cancel that of 

Eq. (|15|) summarises the perturbative jet function of F 2 , namely the IniV enhanced 
contribution to the corresponding coefficient function. It is written in an all-order, 
scheme-invariant way. There are various ways to define the Borel integral in flT5D . In gen- 
eral, different regularizations of the renormalons differ by power terms of the form (|T3|). 
Therefore, once power terms are included it does not matter which particular regulariza- 
tion is chosen ||23|| : difference between different regularizations amounts to redefining Ui. 



Here we choose to define (|T5|) by truncating the corresponding perturbative series 

/^i OO fe + 1 

In J N (Q 2 ; A*!) = ^ £ £ C Kl (\nN) 1 A\ (17) 



A) 



fc=0 1=1 



with A = f3oa s /7r, at the minimal term, i.e. at the point where the terms start increasing 
in magnitude. Even this can be done is several ways, which differ by power corrections. 
One possibility would be to perform the sum over I for a given k and truncate the 
sum over k at the minimal term. We will use another: let us rewrite the sum as in 
resummation with a fixed logarithmic accuracy, 



ex. 



In J N (Q 2 ; /4) = £ 9m{\)A m -\ (18) 
Po m=0 
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where g m (X) are functions of A = A In N, and truncate the sum over m at the minimal 
term. 

It is straightforward to expand the Borel integrand in (fL5l), and derive explicit ex- 
pressions for g m (X). Let us choose \j? F = Q 2 and write g m (\) as a series in 5 = /?i//?o : 



9* 



•(A) * 9m (ty +8g$(\) + Pg${\) + ■■■, (19) 



where first term appears in the large-/?o limit (setting T(u) = 1) and second and on 
are associated with integrating over the two-loop coupling using (|T0|). The expansion in 
powers of S is useful since we will eventually discard higher powers of 5, which enter for 
the first time at rather high orders (note that 5 itself is not small S(Nf = 4) = 0.7392). 
The functions g!$(\) and g$(\) are given by 

9m (A) = r m (\)p m /m\ + 6(m-l)X (p m -a m+1 )/m (20) 
g {2 \\) = 6{m - 1) f m (A)p m _i/ (m - 1)! + 0{m - 2) A (^(m) + lE - l)(p m -i ~ a m ), 

where the coefficients p n (V = A + piA 2 +p 2 A 3 + . . .) are defined through the expansion 
of 

9 

B[P](u) = r(l - u) B[K]{u) = l+p 1 l+p 2 ^ + .... (21) 

In (p0|), aside from a linear term associated with the factorization procedure, the A- 
dependence appears only through the functions 

r m(A) = f A- +2 = r(m - 1) [(1 - A) 1_m - A(m - 1) - 

n=0 l n + 

r m (A) = ^ + lE - 1J r m (A), (22) 

which are definedQ for m > 0. Note that both r m (A) and f m (A) are free of a linear term: 
their expansions start at order A 2 . These functions have two generic properties [pi], |22l [40 



associated with their origin in an integral over the running coupling: they increase 
factorially with m due to infrared renormalons and posses an increasing singularity at 
A = 1, corresponding to Q 2 /N ~ A 2 , where the perturbative treatment looses its validity. 
Terms with 5 2 enter first at NNLL accuracy (m = 2): 

^ (A)= 1 (A + K1-A))' (23) 

The lowest order contribution here is A In J N (Q 2 ; up) = (Cp/Po) | S 2 A 5 In 4 N, where we 
kept only the highest power of In N. Contrary to other terms of the same logarithmic 
accuracy (NNLL) this term has no large numerical coefficients. Normalizing by the 
magnitude of the leading term in In Jn (Q 2 ] ^p) = (C p / ' (3q) \ A In 2 ', and substituting 

10 For m = 0, 1 the limit should be taken carefully, yielding logarithms of (1 — A). 
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AlnN = 1 (for a worst-case estimate) the corrections is of relative size of (5A/2) 2 , i.e. 
~ 0.3% at Q 2 = 10 GeV 2 . In the following we will work with exact NNLL resummation, 
supplemented by estimates of subleading logs (N 3 LL and beyond, up to the minimal term 
in the series). As we will see below, estimates of subleading logs are anyway of limited 
accuracy. So there we will take into account just the first two terms in (|l9l), neglecting 
contributions 0(5 2 ), which are tiny. Note that dealing with N 3 LL and beyond we will 
also be neglecting running-coupling effects associated with four-loop contributions to the 
f3 function, as these corrections are neglected here in the relation between the 't Hooft 
coupling and MS. 

It is useful to relate the Borel representation B[K](u) to the commonly used functions 
A and B. The relation is: 

B[K]{u) — B[A]{u) - uB[B]{u) (24) 

or 

^ A + i^- < 25 > 

If we write B = b x A + b 2 A 2 + b 3 A 3 + . . . with 

9 

~ -7/ ~ II 

B[B](u) = h + b 2 - + b 3 - + . . . . (26) 

then k n = a n+ \ — nb n . Finally, it is straightforward to write the relation with a^ s and 
b^ s defined in the MS renormalization scheme[3 For the first two orders (which are 
known in QCD) we have 

a± = ai = 1 b\ = b\ — —3/4 ko = 1 

a 2 = a 2 = af § ^ _ _ b 2 = b 2 = bf k x = a 2 - h (27) 

a 3 = a 3 - 5a 2 = af + 6f s - 5af s k 2 = a 3 - 2b 2 , 

where 6f* = 0^ //3 3 . The MS coefficients are given by @, % §, [J§, 

^w.A) = | + g-s)|- < 28 > 

- ^ -biM- *) (f " + & ^ " T 

t5 g /Ar 247 1 2 1 // 3 3 A 1 2 \ f-73 5 A \ 

6o s (N c , 8o) = h -7T 2 + — G H — tt ) Cf + I h -Ca C A 

2 v ' M0; 72 6 & VV 32 2 Sd 8 / V 144 2 s 7 

where the purely non-Abelian ingredient in af s is still not known exactly. However, a 
reliable estimate |38| was extracted from the first few moments |37j for N c = 3, allowing 



us to write the above. 



11 Note that our notation for the coefficients of A and B is somewhat different from |3^, |34|. The 
relation is A n = Cp^ n ^ x a^ s and similarly for B n . 
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2.4 Naive non-Abelianization 

As explained in the introduction, our task is to take into account all the perturbative in- 
formation available on In J, namely the large- Nf limit and the first few coefficients known 
in QCD. The first question is, of course, to what extent does "naive non-Abelianization" 
provide an estimate of the coefficients. The answer, as reflected in the numerical values 
of the first few coefficients for N c = 3, 

(29) 

19.5 ±0.2 

Pi 



is quite clear: the terms that are leading in the large-ft limit, are not at all dominant 
for the realistic value ft(iVy = 4) ~ 2. This is also reflected in Table [l] (compare, for 



n 


aZi (large ft) 


C (large ft) 


C s (large ft) 


p n s (large ft) 





+1 





+1 


1 


1 


+1.667 


-0.75 


+2.417 


2.994 


2 


-0.333 


-1.786 


+3.238 


8.006 


3 


-2.071 


+0.425 


-3.347 


22.046 


4 


-1.093 


+7.647 


-31.683 


75.214 


5 


+0.325 


+ 15.591 


-77.629 


343.566 


6 


+0.477 


+ 1.98 


-11.404 


1983.269 


7 


+0.107 


-79.714 


+558.106 


13618.4 


8 


-0.0487 


-223.285 


+1786.23 


107843.3 


9 


-0.0296 


-83.078 


+747.675 


965598.1 


10 


-0.00236 


+ 1319.7 


-13197.2 


9631151.8 


11 


+0.00224 


+4245.12 


-46696.3 


105806143.15 


12 


+0.000694 


+1756.42 


-21076.0 


1268854143.8 



1.667- 
-0.333 + 
-1.786 + 



1.467 



ft ' 
14.714 

ft 
6.6104 



n 


MS 
u n+l 


IMS 
n 


LMS 


MS 





+ 1 





+ 1 


1 


1 


+0.962 


-0.75 


+1.712 


2.289 


2 


+2.232 (2.934) 


+ 1.387 


+0.0115 (0.714) 


4.393 (5.095) 



Table 1: The coefficients of A, B, K, and V defined in (0), (0), (f|) and (|2T|), respectively 
in the large-ft limit (upper table) and in full QCD with Nf = 4 (lower table). In the 
latter, for n = 2 the coefficient in the £ t Hooft scheme is given in parenthesis. 
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example, the n = 2 values in the lower and upper tables). One simply cannot use the 
large /3 limit directly to provide an estimate^ of A, B or /C because these functions do 
not contain renormalons, and thus, there is no reason that the large-/?o terms will have 
large coefficients. In fact, it is known that the series for the cusp anomalous dimension 



A converges in this limit into [|||, [4£ 



sinvrA r(4 + 2 A) 
6r(2 + A) 2 



where the MS scheme is used, and the dots stand for terms that are subleading at 
large (3q. It follows from (12) and (24) that B, is renormalons-free as wellQ. Thus, also 



here "naive non-Abelianization" is not expected to apply. 

The status of V is somewhat different. It does contains renormalons so here the terms 
that are leading in the large-/?o limit are expected to dominate at large orders. This does 
not necessarily imply that "naive non-Abelianization" estimates will be reliable at low 
orders, but it is expected that the general trend of increase characterizing p n in the 
large-/5o limit right from the start, will persist in the full theory. The first few entries in 
the lower Table [I] are consistent with this expectation. 

While the large-/3o coefficients themselves cannot be used to estimate the actual QCD 
ones, more general properties may be deduced. As reflected in Table [I] neither A nor B 
has a regular sign pattern. The sign is flipped every three or four terms owing to the 
sin 7™ function in ( |i~2"|) associated with taking a final-state cut. Another property is that 
in the large-/?o limit the coefficients of A decrease with order while those of B have a 
general tendency to increase. Starting at n ^ 4 there is a clear hierarchy |a n +i| <C \k n \, 
i.e. the physical kernel is dominated by B: k n = a n+ \ — nb n ~ nb n . An even stronger 
hierarchy exists between A and V: |a n +i| ^ \Pn\, as the latter contains renormalons. 

Finally, to make concrete predictions for the subleading logs in (PU|), the knowledge 

of 

Pn ~ a n+ i and p n at higher orders (n > 3) is required. However, from the discussion 
above it is expected that at higher orders |a n+ i| -C \p n \, and thus the knowledge of A to 
higher orders is absolutely inessential. As usual, statements about large orders do not 
necessarily apply for any n > 3. Indeed, as shown in Table |2| the expected hierarchy is 
not realised for n < 3. Yet, the n = 3 entries in Table |l] are reassuring: they suggest 
an order of magnitude difference between p 3 and a 4 . We will therefore orient ourselves 
in the next section to estimating K, and thus V in full QCD while for A we will simply 
use the first known coefficients (a 2 and 63) supplemented by the large-/?o contribution 
(although the latter does not dominate!) for the higher orders. An alternative is to 
neglect a n for n > 4 altogether. 



12 This negative statement does not apply to physical quantities such as the first few moments of the 
structure functions. In this case estimates based on "naive non-Abelianization" are reasonable ^||, 



^3 



3 Although the coefficients of B are known to all orders, a closed-form expression as a function of the 
coupling is not known. 
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n 




Pn 





1 


1 


1 


0.962 


2.289 


2 


2.223 


3.403 



Table 2: The coefficients of B[A]{u) and in QCD with N f = 4. 



2.5 Going beyond the large-/3o limit 

A major advantage of the large-/3o result for K, is that its analytic structure in the Borel 
plane is known. We will use this structure as a basis for modeling the Borel function in 
the full theory. We shall require consistency with the first few orders that are known, 
such that NNLL accuracy is guaranteed. By examining different models we would like 
to check the sensitivity to yet unknown higher-order coefficients, in particular in the 
context of open theoretical problems concerning renormalons. There is no attempt here 
to cover the infinite set of functions that obey the constraints. 

We will assume that the singularity structure of the physical kernel is similar to that 
obtained in the large- (3q limit from (|1^) and (|6|). In particular, in our models there will 
be two renormalons at u = 1 and u = 2, but no higher ones. Thus, non-perturbative 
contributions will be parametrized by fll3"l) . This is not to say that higher renormalon 
singularities do not exist in the full theory, but rather that with the information available 
there is no way we can control their residues. Thus, in our models, similarly to ([T2]). 
B[K](u) does not vanish at u = 1, 2, but it does vanish at all higher integers. 

We will also restrict ourselves to the case where B[lC\{u) is regular at u = 1,2 
so the renormalon singularities of the physical kernel in our models appear as simple 
poles. As mentioned in section 2.2, the nature of the singularities is related with the 
anomalous dimension of higher-twist operators. Here, we keep B[JC](u) renormalon free, 
thus assuming that the anomalous dimension at any twist coincides at large N with that 
of twist two. This issue certainly deserves further investigation in the future. 

We investigate two types of models for B[K](u). The first is based on modifying the 
large-/?o function ( |i~2"D by a multiplicative factor and the second - by an additive term. 
A simple example with a multiplicative factor is 



^r,^n / M 1 SI 11 5, 

£TO| modelo = -— — 
—i a Hi 



i i 

+ 



l-u l-u/2 



x e"^- 



1 + H^u 2 ), 



H {a) = -\a 2 - I - \al + + \at + ^ - \ba 2 - b 2 (31) 

where H^ a ' is fixed to reproduce the exact k\ and k 2 of fl2"7D- Another example (model 
b) is provided by replacing (1 + H^u 2 ) by e H(b) " 2 ; with = H^ a \ Next, a couple of 
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simple examples with additive (regular) terms are 



235 I 



sin nu 



ITU 




+ 



1 



and 



H (c) 



B[JC)(u) 



— f -7T 2 + -af + -5 2 MS 
72 6 2 3 2 2 



M 1 
1 



u/2 
5a 2 - b 2 



+ u [ a 2 - - 



(32) 



model d 



sin ITU [ 1 5„, 

< - es 

2 



ITU 



1-u l-u/2 



_|_ e «(a 2 -f)+« 2 H( d ) _ j 



335 1 



(d) = , ^_2 

72 6 



MS _i_ TMS 

+ 2 9 



^ + 2°3 ' o"2 



-5a 2 - 6 2 



-a, 



o a 2- 



(33) 



The main difference between the two sets of models, is that models c and d keep the 
residues of the renormalon singularities at their large-/3o values, whereas in models a and b 
these residues are modified at NLL by purely non-Abelian contributions, proportional to 
Ca/Po, and, starting from NNLL by correlated multi-gluon emission even in the Abelian 
case. In general the residues should be modified by such contributions, but the way the 
functional form varies is not yet under theoretical control. 

As shown in figure ^ different models extrapolating from the large-/5o limit have very 
different behaviour indeed. The knowledge of the value of B[)C](u) at u = and the first 
two derivative there cannot constrain it away from this point. Note that models with 




Figure 4: Different models for B[K](u). The points describe the calculated function in 
the large-/3 limit fljjD , and the lines describe simple models for this function in QCD, 
inspired by the large-/?o analytic structure, which are consistent with the first three orders 
(^0,1,2) m the expansion of B[K](u). The full lines stand for the following models (from 
bottom to top in the range < u < 1): green - a, blue - b, magenta - c and red - d. The 
dashed line describes model e. Arrows show the points where renormalon poles appear. 



an additive correction (c and d) follow closely the large-/?o functional form away from 
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u = 0, as their residues at u = 1 and at u = 2 are fixed by this limit. On the other hand, 
models with a multiplicative correction (a and b) can have a completely different form 
and even the sign of the first renormalon residue varies. 

In general, the detailed behaviour of B[K](u) away from the origin, has a minor 
significance concerning the value of the physical kernel itself, or J N (Q 2 ; fi 2 F ). To see the 
sensitivity to the behaviour of B[K](u), we computed In J N (Q 2 ; ji 2 F ) with the different 
models described above. The results are shown in Fig. [| order by order at increasing 
logarithmic accuracy. Differences are rather significant. Models a and b tend to have a 
moderate (although not negligible) contributions beyond the NNLL, quite a favourable 
scenario for perturbation theory. On the other hand, models c and d indicate huge 
corrections beyond the NNLL. The resolution of the theoretical questions allowing to 
pick the right model for subleading logarithmic corrections thus becomes crucial for any 
analysis oriented at power accuracy. 




02468 10 2 4 6 

order (N m LL) order (N m LL) 

Figure 5: Order by order calculation of the jet function In J N (Q 2 ; fi 2 F ) for N = 10 and 
for Q 2 = 2 GeV 2 (right box) and Q 2 = 20 GeV 2 (left box). NNLL (m = 0,1,2) are 
exact in all cases. Subleading logs are calculated based on each of the four models (from 
bottom to top): green - a, blue - b, magenta - c and red - d. The minimal term in the 
series is denoted by an 'x' symbol. 

One typical feature of models c and d is that contributions of consequential orders 
are not monotonically decreasing up to the minimal term, but instead, their magnitude 
oscillates. This, however, should not be too surprising given that in the case under 
consideration oscillations in magnitude of consequential terms occur at the first few 
order: the NLL contribution is larger than the LL one, while the NNLL contribution is 
smaller than both. 

As mentioned above, already when constructing a model for B[JC](u) which is con- 
sistent with NLL accuracy in the non-Abelian theory, there are various possibilities. 
Nevertheless, in this case it seems natural to use a multiplicative modification of the 
large-/3o result, such that C = 5/3 at the exponent in ([12]) is effectively replaced for the 
full d2 coefficient (p8|) . In fact, the choice made when promoting the large- (3q result to 
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NLL accuracy is the most significant one. This can be verified by examining the model, 



B[K]{u) 



model e 



sin TTU 



iru 



,«2 U 



+ 



u 



u/2 



+ e u2H(e) - 1 



(34) 



with = H^ a \ where 5/3 is replaced by the full a 2 in the first term, while further 
corrections (NNLL) are incorporated as an additive term. This intermediate step be- 
tween the two classes of models described above, turns out to be quite close to the ones 
with a multiplicative correction, a and h. For example, for N = 10, Q 2 = 20 GeV 2 , upon 
truncation at the minimal term model a yields In J = 0.713, model d yields In J = 0.868 
and model e yields In J = 0.755. 



Indeed, in previous applications of DGE |2T, 22j , the procedure chosen to extrapolate 
from the large-/9o limit and comply with NLL accuracy was the replacement 5/3 — > a 2 . 
This follows the spirit of ]7J who first noted the special role of the "gluon bremsstrahlung" 
effective charge, (defined by the cusp anomalous dimension A) in this context. This nat- 
ural generalization is sufficient at NLL, but not beyond it. In any case, the deeper 
theoretical problem of how purely non-Abelian contributions as well as correlations in 
multiple gluon emission modify the renormalon structure of the Sudakov exponent re- 
mains open. 

With this caution in mind, we will proceed here on good faith with models with a 
multiplicative modification of (0), and in particular, with model a, to perform data 
analysis. A multiplicative modification can always be recast as a modification of the 



coupling which is integrated over by redefining T(u) in (15), and thus it can be regarded 
as the most natural generalization of the large-/?o limit result. Based on Fig. |5| it is also 
an optimistic choice, leading to rather moderate higher-order corrections. 



3 Data analysis 

In the previous section we summarized the state-of-the-art knowledge (and ignorance) on 
the twist-two coefficient function of F 2 at large x. We also have a picture ]T3| concerning 



the structure of higher twist, with a concrete formula for their parametrization (jig). 
We now turn to investigate the implications this has on data analysis. The following 
questions will be emphasized: 

• Are the data and the theoretical description of F 2 consistent? 

• Is there any experimental indication for higher-twist contribution to F 2 1 

• To what extent can one determine the relevant parameters from the data. This 
includes a s , the non-singlet quark distribution ^^(/^f) and the higher twist param- 
eters UJi. 

• What is the significance of Sudakov resummation as compared with the NNLO 
result for the coefficient function? 
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• What is the significance of renormalon resummation in the Sudakov exponent and 
the current uncertainty about subleading logs (beyond NNLL)? 

In order to systematically treat target-mass corrections and simplify the resummation 
procedure, we work here in moment space. For this purpose we compute the Nachtmann 
moments from the available experimental data on F^x^Q 2 ). We also include the con- 
tributions at x — 1 from the proton electric and magnetic form factors. The data we 



use are from a series of SLAC experiments on ep scattering ||50||-||58|| and the BCDMS 
collaboration using \ip scattering |59| . 



Before coming to the fits, let us briefly describe the data and the way we use it. The 
SLAC data cover a region from Q 2 below 1 GeV 2 to 20 GeV 2 and x up to x max defined 
by the minimum value of W 2 = (m p + m^) 2 . The BCDMS data cover a region x < 0.75 
for 20 < Q 2 < 230 GeV 2 . 

The moments are computed in 15 bins of Q 2 from 2 to 105 GeV 2 . There is a region of 
overlap between the SLAC and BCDMS data sets. While small it does allow a check on 
the relative normalisation of the two experiments. The computation involves estimating 
the contribution from the region between the highest x bin at each Q 2 value and x = x max . 
We take a linear extrapolation with an arbitrary uncertainty of 100%. For the SLAC 
data this is a very small correction. The final errors on the values of each moment include 
this uncertainty. 

The physics analysis we apply assumes simply the non-singlet evolution of moments. 
This is equivalent to assuming that high enough moments of F%(x, Q 2 ) are dominated 
by just the u and d valence quarks. This can be checked by taking a standard global fit 
to the DIS data [I] and evaluating the relative contribution to each moment from the 
sea and valence quarks. We find at Q 2 = 10 GeV 2 , for example, that the LO analysis 
of MRST2001 gives valence/total for N = 4 to be over 90% and for N > 5 to be more 
than 95%. We therefore choose to analyse moments N > 5 and keep in mind the valence 
approximation as a source of an overall ^ 5% uncertainty. 

3.1 NNLO based fits 



It is well known [^4J that a NLO analysis of DIS data fails to describe the large-x 
region. Usually the deviations are attributed to higher- twist contributions at low W 2 . 
In the MRST2001 NLO leading-twist analysis, data from W 2 > 12.5 GeV 2 were simply 
dropped. Even including the NNLO corrections which are now available [ 35|] — [ 40[ , though 



helping a little, did not succeed in diminishing the discrepancy at large x. 

This situation is confirmed here by attempting to fit the moments we extracted by a 
NNLO formula, 

F 2 N (Q 2 ) = q N (f4) E N (Q 2 ,n 2 F ) C N (Q 2 ). (35) 
We use here the standard MS factorization and renormalization schemes. The NNLO 
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oc s (M z ) extracted from each moment 



oc s (M 7 ) 



0. 125 — 



0. 105 — 



+ t 



j NNLO only 

■ NNLO -+- pure NNLL re summation 

1 NNLO + full resummation with 
multiplicative form 

X NNLO + full resummation with 
additive form 



+ 



| 



8 9 

N 



10 11 



Figure 6: The extracted value of the coupling from individual moments for Ff in various 
approximations. The symbols represent (from top to bottom): a NNLO fit with no 
Sudakov resummation and three fits with Sudakov resummation, a pure NNLL and two 
renormalon-based models: a - multiplicative modification of the large-/3o limit, and d - 
additive modification. 



evolution factor is 



f 1 f Q2 



« S (Q 2 ) 



cu s (/4) 



[l + r N a s (Q 2 ) + s N a 2 s (Q% 



(36) 

and the NNLO coefficient function with the factorization scale (and renormalization 
scale) set to Q 2 is 

CT L °(Q 2 ) = 1 + c^a s (Q 2 ) + c$a 2 (Q 2 ). (37) 

In (|36D, d N = 7»? /2/?o, r N involves the NLO anomalous dimension 7^ and the NNLO 
anomalous dimension j N . We use the results of Retey and Vermaseren |37j for the 7^ 
for even N and a smooth interpolation of these for odd N. The running of the coupling 
a s (Q 2 ) is approximated by the NNLO (3 function. 

We fit each moment in turn with two parameters, the moment of the valence quark 
distribution ?jv(a*f) an d af s (Mg). The resulting values of the latter are shown in Fig. ^| 
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NNLO only 





Figure 7: Experimental moments N = A through 11 of with two NNLO descriptions 
for each N: one (dashed line) based on fitting ^(/^f) an d af s (Mf ) and the second (full 
line) based on extrapolating the standard MRST fit to low W 2 . 

(upper box) showing two features: (i) large value of the coupling and (ii) a marked 
increase as N increases. The MRST2001 NNLO fit with W 2 > 12.5 GeV 2 prefers a much 
lower value, af s (M|) = 0.1155. In Fig. ^ we show two NNLO descriptions, one with the 
pairs of parameters just described and the other corresponding to the extrapolation of 
the MRST2001 NNLO description to low W 2 . We conclude, with no surprise, that the 
leading- twist NNLO formula cannot satisfactorily describe the large N moments of F%- 

3.2 Resummation based fits 

Next consider the additional effect of Sudakov resummation. This is done by modifying 
eq. ([J7|) to include higher-order terms which are enhanced by powers of IniV. We use, 
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NNLO ■+- Resummation model (a), A. =0.4, N - =5 



uncertainty 
► extrapolation 





lO — 




Figure 8: Experimental moments N = 4 through 11 of and a perturbative description 
based on Sudakov and renormalon (model a) resummed coefficient function, matched 
into the NNLO result. Here a^(Mf) = 0.1135 (this is the best fit value). The full line 
represents the region on which the fit is based whereas the dashed line represents an 
extrapolation to lower W 2 . Dotted lines represent the resummation ambiguity based on 
(±) the size of the minimal term. 



the so-called log-R matching, namely 

\nC N (Q 2 ) = In C N N Nh °(Q 2 ) + In J N (Q 2 ) - \nC°^ p (38) 

where In Jn{Q 2 ) is given by (O) or ( |T8"D with the factorization scale set to Q 2 . The sub- 
tracted term accounts for log-enhanced contributions up to order a 2 which are included 
in both the NNLO coefficient function and in the Sudakov resummed one, 



^overlap = CF 



h 2 + { lE -b 1 )L^A (39) 



+ ( ^ + ( + ^(a 2 - h) ) L 2 + ( lE (a 2 -b^-h + ^ + ^j Lj A 2 
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Figure 9: Values of the parameters (In^f) f° r = 2,5,10 and 20 GeV 2 . The differences 
between the values for the two cases (i) W 2 < 12.5 GeV 2 excluded with NNLO evolution 
only and (ii) low W 2 data included with NNLO + Sudakov resummation indicate the 
modifications required to the large N moments due to resummation effects. 



We apply the two following criteria in selecting data for the fits. Firstly, we fit 
moments iV > 5 , to ensure non-singlet dominance. Secondly, we limit A = AlnN = 
(j3 a s /7i)hiN < A max , to ensure the validity of the perturbative treatment. We choose 
Amax = 0.4. Since A depends on one of our fit parameters, cn s (M§), the number of data 
points included in a fit could vary with this parameter. To avoid this we choose the data 
to be fitted for A max evaluated for the MRST2001 NNLO default value, a s {M 2 z ) = 0.1155. 
The moments of the experimental data going into the fits are given in Table Note 
that we chose this same data set for the pure NNLO fits for the sake of comparison. 

As in the NNLO fit we have two free parameters, qN(p F ) and af s (M§). We choose 
three variants for Sudakov resummation. The first is the standard, NNLL accuracy 
resummation, and the other two include renormalon resummation according to models 
a and d discussed in section 2.5. Contrary to fixed logarithmic accuracy, renormalon 
resummation is oriented at power accuracy. The renormalons are regularized here simply 
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iV 


/-\'2 
*Vmin 


no. of points 


5 


2 


15 


6,7,8 


2.8 


14 


9,10 


3.6 


13 


11 


4.5 


12 



Table 3: The experimental data points included in the fits. For all N values Q^ax = 105 
GeV 2 . The number of data points with different Q 2 values is indicated in the last column. 

by truncating the series (|I3) at the minimal term. 

The results are shown in Fig. |6|. In contrast to the NNLO case, the values of the 
coupling are practically independent of N. This is true for any of the renormalon-based 
models as well as for pure NNLL resummation. Thus, Sudakov resummation proves to 
be important: it provides the ingredient which is missing in the fixed-order perturbative 
description discussed above. 

However, as shown in Fig. the extracted value of a s (M§) varies significantly de- 
pending on the chosen model. In pure NNLL resummation, typically a s (M|) = 0.1155 ± 
0.00100 while with renormalon resummation at a (M§) varies between 0.1135 ± 0.0010 
(model a) and 0.1070 ± 0.0010 (model d). Since in any case NNLL accuracy is guar- 
anteed, these differences are due to sub-leading logs, N 3 LL and beyond. Note that the 
result quoted for model a actually represents a large class of models with a multiplicative 
modification of the large- (5q all-order result. The details of this factor have a minor effect 
of the result. Taking the optimistic approach that these multiplicative models represent 
well the structure of the Sudakov exponent, the central value of the strong coupling 
from this analysis is 0.1135. Clearly, more theoretical work is required to reduce the 
uncertainty. 

Fig. |S] shows the results of fitting q N based on model a for a fixed value of a s (M|) = 
0.1135. The plotted curves are almost identical to those obtained in a similar fit where a s 
is free, since the variations of the best-fit value of a s as a function of N are small (Fig. |6|). 
Note that in summing the series (|18"D, the asymptotic expansion is terminated at the 
minimal term. A sensible estimate of the corresponding power-suppressed ambiguity is 
simply the size of the last term included in the sum. Fig. [8] shows the magnitudes of this 
ambiguity together with extrapolations of the fits to regions excluded from the fitting 
procedure. Under the assumption taken here that the dominant higher-twist terms are 
those associated with mixing with the leading twist, this power-suppressed ambiguity 
should provide an estimate of the magnitude of the higher twist. As clearly seen in 
Fig. |H this is a small correction even in the highest moment. These results imply that 
the low W 2 data, i.e. W 2 < 12.5GeV 2 can be successfully described upon improving the 

14 The error quoted, like the bars in Fig. |[ reflects just the propagated experimental error (taking 
into account the uncertainty in extrapolating the data to large x). It does not include any estimate of 
theoretical uncertainty, for example due to factorization or renormalization scale dependence. 
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standard NNLO description by including Sudakov resummation. The resulting values of 
a Ms (_/y/|) are stable with N and are close to the value extracted from the MRST NNLO 
analysis where a cut on low W 2 is imposed. 

The other parameter ?at(a*|') is shown in Fig. § for four values of /i^ for comparison 
with the corresponding values expected from the MRST NNLO analysis. For the most 
part the ranges of x being explored in this analysis of the high moments correspond to 
ranges of W 2 below the value 12.5 GeV 2 excluded by the MRST analysis. Thus Fig. |9] 
is an indication of how much the large x parton distributions would have to be modified 
to take account of this resummation effect. 

3.3 Fitting the higher twist 



NNLO -+- Mult. Resum. -+- Higher Twist 




0.05 0. 1 

C HT (GeV 2 ) 



Figure 10: The total \ 2 as a measure of the quality of the fits (upper box) to individual 
moments when C HT is varied and (^(/-t^) and a s (M|) are free parameters, and the 
corresponding best-fit value of a s (M§) (lower box). Model a is assumed. 

Finally, according to eq. (|13|), we include the simplest model for higher- twist correc- 
tions through an extra factor J NF (NA 2 /Q 2 ) where 

J-(iVA7Q 2 )=expj^j, (40) 
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with C HT = ~ UJ ^J^ Here we neglect subleading corrections proportional to u^. 

We perform fits on the same set of moments as before, but with the extra parameter 
C HT characterising the additional Q 2 /N dependence from the higher twist. Naturally, 
there is a compensation between the magnitude of the parameters a: s (M§) and C HT which 
is shown in Fig. |10| in the lower box. 

For a wide range of C HT there is little difference in the quality of the fits, see the 



upper box in Fig. [10]. The conclusion is clear: there is no phenomenological evidence for 
a non-vanishing twist-four (or higher twist) contributions going beyond the resummed 
perturbation theory. On the other hand a positive C HT would decreases the already 
'lowish' value of a.. 



4 Conclusions 

In this paper we investigated the large-x limit of the structure function F 2 (x, Q 2 ), concen- 
trating on the interplay between Sudakov resummation, renormalons and higher twist. 
The main assumption we make [T3[ is that the dominant higher-twist contribution in this 



region is related to the mixing of higher-twist operators with the leading twist. Con- 
sequently, our investigation focuses on constraining the twist-two coefficient function at 
large-x. 

Our first point is that the physical kernel dlnF 2 N jd\wQ 2 receives log-enhanced con- 
tributions from a single source, the jet function depending on W 2 = Q 2 (l — x)/x. This 
property is a consequence of factorization and it distinguishes structure functions from 
less inclusive quantities, such are event-shapes, which have an additional observable- 
dependent sensitivity to large angle soft emission. 

Much is known about the structure of the Sudakov exponent of F 2 from general 
considerations, from the large-/?o limit and from fixed-order calculations in dimensional 
regularization. Nevertheless, the investigation of Sec. 2.5 reveals that subleading logs, 
N 3 LL and beyond, which are not yet under theoretical control, may have a large ef- 
fect. Further theoretical input on the all-order structure of the exponent is required to 
constrain these contributions. 

While a pure NNLO formula without Sudakov resummation does not describe the 
large-x data, once Sudakov resummation (NNLL) is employed the purely perturbative 
description of F 2 is fully consistent with the data. This is reflected in the good fit of 
the Q 2 dependence of each of the moments and in the stability of the extracted a s 
as a function of N. This conclusion holds independently of whether renormalons are 
resummed and of what one assumes about subleading logs. 

The theoretical uncertainty concerning the renormalon structure of the Sudakov ex- 
ponent beyond the large-/?o limit, as reflected in the models investigated in Sec. 2.5 
and 3.2, translates into an uncertainty of order ±6% in the extracted value of the cou- 
pling. Thus, taking a conservative approach one would assign a rather large uncertainty 
to a s from scaling violation of F2 (Q 2 ) with N <: 5. A more optimistic approach would 
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be to assume that a multiplicative modification of the large- (5q result, such as model a 
in Eq. ([H]), represents well the all-order structure of the exponent. In this case the total 
effect of subleading logs, N 3 LL and beyond, is similar to the NNLL contribution (see 
Fig. H). As a result the central value of ct s (M|) changes from the pure NNLL result of 
0.1155 to 0.1135. 

A priori, higher twist is expected to be important for F 2 at large x. However, as it 
stands, there is no experimental evidence for a non-vanishing higher-twist contribution 
going beyond the resummed perturbation theory. This is consistent with (at least, it does 
not contradict) the assumption of [13|] that the dominant higher twist at large-x is the 
one related to mixing with the leading twist. On the other hand, in the experimentally 
accessible range and excluding the very small Q 2 /N region (Q 2 /N ~ A 2 ), the natural 
size of the higher twist under the above assumption is still significantly smaller than the 
theoretical uncertainty. 

In addition to solving the theoretical problems mentioned above, progress in this field 
requires additional data. Although the SLAC and BCDMS data are reasonably precise 
it is desirable to extend the range they cover. Ideally, to verify or falsify Eq. ([!]) one 
would need data covering a large range of N and Q 2 with a fixed Q 2 /N. 
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